Probing CPT violation in neutrino oscillation: A three flavor analysis 
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We have studied CPT violation in neutrino oscillation considering three flavor framework with mat- 
ter effect. We have constructed a new way to find the oscillation probability incorporating CPT 
violating terms without any approximation. Then CPT violation with atmospheric neutrinos for a 
magnetized iron calorimeter detector considering the muons (directly measurable with high resolu- 
tion) of the charge current events has been studied for zero and nonzero S13 values. It is found that 
a potential bound of Sbs2 ^ 6 x 10^^* GeV at 99% CL can be obtained with 1 Mton.year exposure of 
this detector; and unlike neutrino beam experiments, there is no possibility to generate 'fake' CPT 
violation due to matter effect with atmospheric neutrinos. The advantages of atmospheric neutrinos 
to discriminate CPT violation from CP violation and nonstandard interactions are also discussed. 

PACS numbers: 14.60.Pq 



I. INTRODUCTION 

The invariance of the product of charge conjugation 
(C), parity (P) and time reversal (T) CPT is one of the 
most fundamental symmetries in physics and it is indeed 
related intimately with Lorentz invariance. The CPT 
theorem [l| states that any local quantum field theory 
(QFT) which is Lorentz invariant and has a Hermitian 
Hamiltonian must have CPT symmetry. The general 
theoretical proof of CPT invariance in particle physics 
along with accurate experimental tests of CPT violation 
(CPTV) provides an attractive candidate signature for 
non particle physics such as string theory [3, y] . In par- 
ticular, the assumptions needed to prove CPT theorem 
are invalid for strings which are extended objects. More- 
over, since the critical string dimensionality is greater 
than four, it is possible that higher dimensional Lorentz 
invariance breaking would be incorporated in a realistic 
model. O.W. Greenberg has shown that CPT violation 
implies violation of Lorentz invariance, but, CPT invari- 
ance is not sufficient for out-of-cone Lorentz invariance 



Consequences of CP, T and CPT violation in neu- 
trino oscillation have been discussed in |21|. Briefly, for 
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j/q — > vp flavor oscillation probabilities Pap at a distance 
L from the source is the following. If 



then CP is not conserved. If 



/3^a, 



then T-invariance is violated. If 



(1) 



(2) 



or 



Pc,p{L) ^ PpdL), 13^ a, (3) 

Po.a{L) ^ Paa(i), (4) 



then CPT is violated. Also it has been shown that mat- 
ter effects give rise to apparent CP and CPT violation 
even if the mass matrix is CP conserving. 

At the experimentally accessible energies, the signals 
for Lorentz and CPTV have been described by a class of 

nrt 

extensions of the Standard Model (SME) [5-8]. The low 
energy probes of these new physics have been explored 
in many experiments with current technologies [9|. A 
summary on measured and derived values of the coeffi- 
cients for Lorentz and CPTV in SMEs have been nicely 
tabulated in [10| for matter, photon and gravity sectors. 
The SME framework predicts several unusual phenom- 
ena in neutrino oscillation, among which are neutrino- 
antineutrino oscillation, directional dependence, uncon- 
ventional energy behavior ll|. These effects have been 
explored with recent experimental data at MINOS [18 1 
and at LSND 19| and new bounds have been obtained. 



The effect of perturbative Lorentz and CPTV on neu- 
trino oscillation has also been studied in 12| with leading 
order corrections arising from renormalizable operators 
including the above phenomena. More consequences of 
different SMEs of Lorentz and CPTV in neutrino oscil- 
lation can be found in 



13 



15J for atmospheric solar and 
baseline experiments as well as in [16|, |l7| for ultra-high 
energy neutrino experiments. 

With an effective-theory approach, a plausible spon- 
taneous CPT violating minimal extension of standard 
model that are also Lorentz invariance violating (LV) has 
been done in [6| for Dirac particles and in [7| for Majo- 
rana particles. Our analysis applies for both Dirac or 
Majorana neutrinos. 

In these theories the Lagrangian for a fermion to the 
lowest order in the high scale can be written as 

C = i'09^7^V - ™V''0 - ^m'^7^V' - SpV'757'''0 , (5) 

where, A^j^ and B^ are real numbers. The terms contain- 
ing Aij, and i?^ are clearly Lorentz invariance violating. 
The effective contribution from these terms to the neu- 
trino Lagrangian can be parametrized as 



/:r^^ = ^£&f7''-f; 



(6) 



where, h^ are four Hermitian 3x3 matrices corresponding 
to the four Dirac indices /x and a, /3 are flavor indices. 
The effective Hamiltonian for ultra-relativistic neutrinos 
with a definite momentum p is then 



H; 



2p 



+ h 



(7) 



where M is the neutrino mass matrix and b = 5o in the 
CPT conserving limit. If we choose a preferred frame 
in which the Cosmic Microwave Background Radiation 
(CMBR) is isotropic, then the rotational invariance im- 
plies no directional dependence for b [7|. 

In presence of LV, one can obtain the same effective 
Hamiltonian by considering a modified dispersion rela- 
tion for neutrinos, E^ = F{p,m). Using the rotational 
invariance in the CMBR frame and demanding Lorentz 
invariance at low energy, this dispersion relation can be 
written as [20| 



/ 



(3) 






(8) 



where, /'■"^'s are dimensionless quantities. Epi is the 
Planck scale of energy where the Lorentz invariance is ex- 
pected to be broken. For ultra relativistic neutrinos with 
fixed momentum, the above dispersion relation becomes 



E = p- 



m 



b--- 



(9) 



with b = Epif^^^ /2 as the leading CPT violating con- 
tribution. For three flavors, it leads to the same effective 
Hamiltonian as in Eq. [71 

The CPTV in neutrino oscillation was first proposed 
for two flavor case in [7| . The typical frequency of neu- 
trino oscillation is Am'^/{2E). For atmospheric and long 
baseline neutrinos, it can be as small as 10^^^ GeV. If 
the accuracy of the oscillation frequency is 10%, then one 
can naively estimate the CPTV parameter to the order 
of 10-23 GeV. 

The CPTV in two flavor formalism has been studied 
for the future atmospheric neutrino experiment at a mag- 
netized iron calorimeter (ICAL) detector 22| and for long 
baseline (735 km) experiment with a typical neutrino fac- 
tory 23j . In 23| , it has been shown that the interference 
between the CPT violating interaction and CPT— even 
mass terms in the Lagrangian can lead to a resonant en- 
hancement of the oscillation amplitude and this may lead 
to 'fake' CPT violation. In [24], it has been shown that 
for hierarchical mass spectrum the upper bound of the 
neutrino antineutrino mass difference can be achieved in 
a neutrino factory is Imp — rh^l ^ 1.9 x lO""*^ eV. 

The LSND result 12a] when combined with the solar 
and atmospheric neutrino observations indicated three 
distinct neutrino mass squared differences. Then it was 
proposed that the CPT violating effects may be large 
enough to make the neutrino and antineutrino spectra 
significantly different 26|, |27|. However, this fact was 



found not to be viable when combined with other neu- 
trino experiments 



28j , and the subsequent data of oscil- 



lations corresponding to AttIq in antineutrinos at Kam- 
LAND [29| ruled it out. The CPT violation is not re- 
quired to explain any neutrino oscillation data if the 
LSND results are ignored in the light of the negative re- 
sults of MiniBooNE [30| that explore the same parameter 
space. However, the current uncertainties in the mea- 
surements of Am2i and Am^2 Ml) allow the possibility 
to study CPT violating effects in neutrino oscillations. 



which may be observed or constrained at the future high 
precision neutrino oscillation experiments. 

From the present experimental data, in Ref. |3l| . it is 
shown that Sb ^ 1.6 x 10~^^ GeV with solar and Kam- 
LAND data, and, in Ref. [33, (5& < 5 x lO^^s GeV with 



full atmospheric and K2K data for two flavor analysis. 

The three flavor analysis of CPTV is mainly needed 
to account properly the matter effect and CP violating 
phases. If b is not diagonal in mass basis, one needs to 
introduce an unitary matrix Ub to diagonalize it. The 
matrix Ufc then needs three angles ^612, ^h23; ^bi3 ^-i^d 
six phases for a complete parametrization and it can be 
written as 

Vb{{Obij}; {0bi}; {ati}; St) 
= diag(l, e^*-, e'^"^) • UcKAf ({^fc^ }; Sb)- 



diag(e 



IQbl 4062 plCtba 



,e 



,e 



') 



(10) 



where, abi, ab2, ctba are the Majorana phases and will 
not have any contribution. To diagonalize the effective 
Hamiltonian, the mixing matrix appears in term of a to- 
tal of six mixing angles (6112, 6*23, ^13, ^6127^623- ^613) and 
four phases {Scp,Sb,<pb2,(f'b3)- This makes the marginal- 
ization of the data practically impossible with normal 
CPUs. 

In [33[, the authors have analyzed CPTV for three 
flavor case and have treated the effect of the CPT vio- 
lating term as a perturbation parametrized by a dimen- 
sionless auxiliary parameter e = 0.1. Finally, they iden- 
tify the combinations of CPT violating parameters that 
contribute to the probabilities to leading order in e and 
compare the signals in different channels to estimate the 
extent to which these CPT violating combinations can be 
constrained or identified in future long baseline experi- 
ments. 

In this paper, we assume b matrix to be diagonal 
in mass basis. Then, there arises only two parameters 
6b2i = b2 — bi and (5632 = 63 — &2- we have introduced 



a new method following Ref. 3j|, where diagonal b ma- 
trix directly comes into the picture. So, in this formalism, 
there is no need to transfer b to flavor basis. Here, we 
have also assumed a preferred frame in which the CMBR 
is isotropic, then the rotational invariance implies no di- 
rectional dependence for b [7|. 

Then with this three flavor formalism, we have stud- 
ied CPTV with atmospheric neutrinos for a magnetized 



ICAL detector proposed at the India-based Neutrino 
Observatory (INO) [35] considering the muons (directly 
measurable quantities at ICAL). In this type of detector, 
one can separate neutrinos and antineutrinos due to the 
magnetic field. 



II. NEUTRINO OSCILLATION IN MATTER 
WITH CPT VIOLATION 

We have incorporated CPTV interactions in a new 
way where two CPTV terms come directly in the oscilla- 
tion probability without any approximation. We derived 
this formalism from the original work in 3J| . We discuss 
it here for any number of generations. 

The neutrino flavor eigenstates |i'a(a = e,/i, r...) and 
mass eigenstates \i'i{i = 1,2,3...) at time t = are con- 
nected by a unitary transformation, 



Wa) = y^^UaiWi)- 



(11) 



Then for a relativistic neutrino beam with energy E, 
the standard amplitude A and probability P for v^ ^ ^^ 
transition after time t in vacuum are 

Aii^c. ^ f^/j) = X! Ua,^exp{~-im1t/E)U}p, (12) 



P{va -^ yp) = l^(i^a -^ i^/j)!^; 



(13) 



where, rrii are the mass eigen values and since neutrinos 
are ultra relativistic t/E = L/E in units h = c = 1. 

To treat neutrino oscillation in matter, let us introduce 
an arbitrary state vector in neutrino flavor space. 



E^^WI'^') 



(14) 



Now, at time i = 0, consider the state as I'a- Then 
'0i(O) — Uai and the transition amplitude is 



Aiiy^^iyf3) = Y,UlMt)- 



(15) 



Then time evolution in presence of matter and CPTV 
term becomes 



id^j{t)/dt = [m'j/{2E) + bj]ijj{t) 

+ Y.V2GN,UejUlMt) 

k 

= Hkjipkit) 



A{Va -^Vp)=^^ UaiXijUJp 



(16) 



assuming b to be diagonal in mass basis. Here, N^. and 
G are the electron number density of the medium and 
Fermi constant, respectively. This equation is for neutri- 
nos and for antineutrinos, the sign of y/2GNe and bi will 
be reversed and U will be replaced by U* . 

The problem of propagation is therefore to diagonalize 
the matrix H defined in Eq. 1161 It can be diagonalized by 
defining a new basis states \v[) = V^J\Va) = ^iLc^Qjl^j)i 
where V is the unitary matrix. If Mf /{2E) are the set of 
eigen values of H , then the solution for a uniform medium 
is given by replacing rrii by Mi and U by V ^ respectively. 

However, one can find the solution without explicitly 
finding V in the following way. For n generation, Eq.. 
1161 has n independent solution for the row vectors ^j (t) . 
At time i = 0, we choose the set of solutions ^/j^' {i — 
1, • • • n) for the row vectors. These are pure mass eigen 
states. 



^f\t^0)^5,, 



(17) 



When these row vectors are assembled into a. n x n 
matrix X according to 



X,,(i)-^f(i), 
then X follows the matrix equation 



(18) 



idX/dt = XH 

with boundary condition X(t — 0) 
an analytical solution is possible. 



(19) 
1. If A^e is constant, 

(20) 



X{t) = exp(-iHt). 

The row i oiX matrix describes the state \i'i) which is a 
mass eigen state at the starting point (t — 0) and column 
j describes the amplitude for evolving into the mass eigen 
state \i'j) at time t. Then the transition amplitude in 
matter is 



where, 



X 



E 

k 



H 



j#fe- 



{2EH - Mfl) 



exp 



' 2E 



(21) 



(22) 



with SM^^ = Ml - Mf. If one subtracts m\/{2E) and 
&i, which will not alter the oscillation probability, there 
will appear Am^^ and Sbji = bj — bi in the final formula. 
We have developed a numerical program considering 
the above formalism using Preliminary Reference Earth 
Model (PREM) [36| for the density profile of the Earth 
to find X and then the oscillation probabilities. We have 
checked that when Sbji = 0, the result matches exactly 
with the other methods of finding oscillation probabilities 
in matter. 



III. 



CPT VIOLATION WITH ATMOSPHERIC 
NEUTRINOS 

A. The chi-square analysis 



We have studied the atmospheric neutrino data for 
the magnetized ICAL detector considering the muon en- 
ergy and direction (directly measurable quantities) of 
the events in context of CPT violation. Due to rela- 
tively heavy mass of muon, it looses energy mostly via 
ionization and atomic excitation during its propagation 
through a medium. Since ICAL is a tracking detector, it 
gives a clean track in the detector. The muon energy can 
be measured from the bending of the track in magnetic 
field or from track length in case of fully contained event. 
The direction can be measured from the tangent of the 
track at the vertex. From the GEANT 



37| simulation 



of ICAL detector, it is clear that the energy and angular 
resolutions of the muons are very high (4-10% for energy 
and 4-12% for zenith angle) and negligible compared to 
the resolutions obtained from the kinematics of the scat- 
tering processes. In our analysis, we neglect the effect 
arisen due to the properties of the detector. 

A new method for migration from true neutrino en- 
ergy and zenith angle to muon energy and zenith angle 



has been introduced in [38j and subsequently used in 



39- 



43|. Here we used NUANCE-v3 '44*] for generating the 



events. The addition of the hadron energy to the muon 
energy, which might improve the reconstructed neutrino 
energy resolution, is not considered here for conserva- 
tive estimation of the sensitivity. It would be realistic in 
case of GEANT-based studies since the number of hits 
produced by the hadron shower strongly depends on the 
thickness of iron layers. 



The x^ is calculated according to the Poisson proba- 
bility distribution. The binning the data is made in 2- 
dimensional grids in the plane of logj^Q E - L°-^. For each 
set of oscillation parameters, we integrate the oscillated 
atmospheric neutrino flux folding the cross section, the 
exposure time, the target mass, the efficiency and the two 
dimensional energy-angle correlated resolution functions 
to obtain the predicted data for the x^ analysis. We use 
the charge current cross section of Nuance- v3 144| and the 
Honda flux in 3-dimensional scheme [45|. This method 
has been introduced in [SSl , but the number of bins and 



resolution functions have been optimized later 391-142 [. 
The method for migration of number of events from neu- 
trino to muon energy and zenith angle bins, the number 
of bins, the systematic uncertainties, and the cuts at the 
near horizons are described in 39|. Both theoretical and 
experimental data for x^ analysis have been generated 
in the same way by migrating number of events from 
neutrino to muon energy and zenith angle bins using the 
resolution functions, which has been used in our previous 
work 1411 ■ 



We marginalize the x^ over the oscillation parameters 
Am§2, ^23, ^13, ^CPi '5^21 and 6b32 for both normal hi- 
erarchy and inverted hierarchy with i^s and Ds separately 
for a given set of input data. Then we find the total x^ 



as x^ = xl 



/Cu ■ 



We have chosen the range of Arn^ 



32 



2.0 - 3.0 X 

lO-^eV^, 6I23 = 37° - 54°, e^ = 0° - 12.5° and ,5621,32 = 
- 5 X 10-23. Am?.^ = 7.06 - 8.34 x IQ-^eV^ and 
9i2 = 30.5° — 40°. However, the effect of Am^i comes in 
the subleading order in the oscillation probability when 
E ^ GeV and it is marginal. We have set the inputs of 

^32 



Am§2 = -2.5 x IQ-^cV^, 6I23 = 45° and 6cp = 0. 



B. Result and discussions 

We have studied the atmospheric neutrino oscillation 
for 1 Mton.year exposure (which is 10 years run of 100 
kTon) of ICAL. The bounds on Sb2i and (5632 are pre- 
sented in Fig. [H It is found that the bounds are stronger 
for 013 — than its nonzero values. This can be un- 
derstood in the following way. The atmospheric neutri- 
nos covers wide ranges of E and L. The matter effect 
only plays role for neutrino with normal hierarchy and 
for antineutrinos with inverted hierarchy. The matter 
resonance occurs in a limited zone oi E — L plane when 
013 is nonzero. The resonance zones squeeze rapidly with 
decrease in ^13 values, while the CPTV contribution is 
independent of energy, baseline and ^13 values. So, the 
CPTV effect is only smeared out to some extent due 
to matter effect and bounds are weaker. It is very im- 
portant to note that no 'fake' CPTV signal is arisen 
with atmospheric neutrinos, which may be the possibil- 
ity in case of baseline experiments. We have shown that 
a strong bound of ,5632 < 6 x IQ-^^ GeV at 99% CL 
can be obtained at a magnetized ICAL detector at INO, 
which is much stronger, almost one order of magnitude 
higher than those estimated till now for future experi- 
ments. However, the bound on Sb2i is very poor, which 
is expected to improve in an analysis with solar neutrino 
data. 

It can be understood that CP and CPT violation ef- 
fect can be very cleanly discriminated in atmospheric 
neutrino oscillation data. The effect of ^13 and 5cp ap- 
pears dominantly neither in atmospheric neutrino oscil- 
lation nor in solar neutrino oscillation, but appears as 
subleading in both cases. These are observable in range 
of E ^ 1 GeV for atmospheric neutrino, where solar and 
atmospheric neutrino oscillations couple. See Fig. 1 and 
2 in Ref. [40[. It has also been shown in Fig. 7 in Ref. 



40| that ii E,y > 2 GeV, the sensitivity to CP phase 
is very negligible. The sensitivity to CPTV comes from 
the sensitivity of Am^2 of the experiment. From our pre- 



vious work [42], it can be understood that sensitivity to 
Arn^2 '^ill ^^ot be lost significantly if one consider events 
with Ey > 2 GeV. So, E^ — 2 GeV can be considered 
roughly the boundary line for CP and CPT violation 
studies with atmospheric neutrinos. 

The nonstandard interactions may mimic as a signa- 





5 b32 (GeV) 



5 b2, (GeV) 



FIG. 1: The bounds on 5632 (left) and ^621 (right) for input 6*13 = 0° and 7.5°. The curves for 5621 with e^is = 0° and 7.5° have 
overlapped on each other. We have set other oscillation parameters at their best-fit values and considered inverted hierarchy as input. 
The marginalization has been carried out over whole allowed ranges of all oscillation parameters (discussed in the text). 



ture of CPT violation. Hovirever, this can also be sep- 
arated out to an extent using atmospheric neutrino os- 
cillation at a magnetized detector. In case of nonstan- 
dard interactions, its contribution comes through mat- 
ter resonance which occurs for some particular zones in 
L — E plane of neutrinos with normal hierarchy and for 
antineutrinos with inverted hierarchy. These zones does 
not move significantly with the change of oscillation pa- 
rameters for their present ranges of the uncertainties. On 
the other hand, CPTV effect comes over whole L — E 
plane. So, again separating the resonance zones, one can 
separate out the CPT violation from nonstandard inter- 
actions. 

The separation of CP violation, CPT violation and 
nonstandard interactions can be done for atmospheric 
neutrinos since it covers wide ranges of E and L and 
contains both neutrinos and antineutrinos. This is a ma- 
jor advantage of atmospheric neutrino experiment with 
a magnetized detector over neutrino beams. 

IV. CONCLUSION 

We have studied CPT violation in neutrino oscillation 
in a full three flavor framework with matter effect con- 



structing a new way to find the oscillation probability 
so that CPT violating terms are incorporated in the os- 
cillation probability formula without any approximation. 
We carried out this study with atmospheric neutrinos for 
a magnetized iron calorimeter detector considering the 
muons (directly measurable with high resolution) of the 
charge current events. We have shown that a stringent 
bound of (5632 ^ 6 X 10^24 Qgy ^^ 99% ,j^ ^^^ ^3^ ^Y}_ 

tained at a magnetized ICAL detector with 1 Mton.year 
exposure, which is much stronger, almost one order of 
magnitude higher than those estimated till now for fu- 
ture detectors. The advantages of atmospheric neutrinos 
to discriminate the effect of CPT violation, CP violation 
and nonstandard interactions are also discussed. 
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